UT-STPD-1/02 
FISIST/05-2002/CFIF 



Yukawa Quasi-Unification 

M.E. G6mez(i)*, G. Lazarides(2)t and C. Pallis^^^ 

^-^^ Centra de Fisica das Interacgoes Fundamentais ( CFIF), Departamento de Fisica, 
Instituta Superior Tecnico, Av. Rovisca Pais, 1049-001 Lisboa, Portugal. 
^"^^ Physics Division, School of Technology, Aristotle University of Thessalaniki, 

54124 Thessalaniki, Greece. 
Scuala Internazianale Superiare Di Studi Avanzati (SISSA), 
Via Beirut 2-4 , 34013 Trieste, Italy. 



Abstract 



We construct concrete super symmetric grand unified theories based on 
the Pati-Salam gauge group 5C/(4)c x SU{2)l x SU{2)ji which naturally 
lead to a moderate violation of 'asymptotic' Yukawa unification and thus 
can allow an acceptable 6-quark mass even with universal boundary con- 
ditions. We consider the constrained minimal supersymmetric standard 
model which emerges from one of these theories with a deviation from 
Yukawa unification which is adequate for /j, > 0. We show that this model 
possesses a wide and natural range of parameters which is consistent with 
the data on b 57, the muon anomalous magnetic moment, the cold 
dark matter abundance in the universe, and the Higgs boson masses. The 
lightest supersymmetric particle can be as light as about 107 GeV. 
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I. INTRODUCTION 



The most restrictive version of the minimal supersymmetric standard model (MSSM) 
with gauge coupling unification is based on radiative electroweak breaking with universal 
boundary conditions from gravity-mediated soft supersymmetry (SUSY) breaking and 
is known as constrained MSSM (CMSSM). It is desirable to further restrict this model 
by assuming that the t-quark, 6-quark and r-lepton Yukawa couplings unify 'asymptoti- 
cally', i.e., at the SUSY grand unified theory (GUT) scale Mqut ~ 2 x 10^'^ GeV. This 
assumption (Yukawa unification) naturally restricts the t-quark mass to large values 
compatible with the data. Also, the emerging model is highly predictive. 

Yukawa unification can be achieved by embedding the MSSM in a SUSY GUT with 
a gauge group containing SU{4:)c and SU{2)r. Indeed, assuming that the electroweak 
Higgs superfields hl^, h'^ and the third family right handed quark superfields f^, h'^ 
form SU{2)r doublets, we obtain the 'asymptotic' Yukawa coupling relation ht = hh 
and, hence, large tan/3 ~ mt/rrib- Moreover, if the third generation quark and lepton 
SU{2)l doublets (singlets) gs and /s {W and r^) form a SU{A)c 4-plet (4-plet) and the 
Higgs doublet h']^ which couples to them is a S'f/(4)c singlet, we get = hr and the 
'asymptotic' relation = follows. The simplest GUT gauge group which contains 
both SU{A)c and SU{2)r is the Pati-Salam (PS) group Gps = SU{A), x SU{2)l x SU{2)r. 
Higher groups with the same property are 5*0 (10) or Eq. 

One problem, which is faced in trying to incorporate Yukawa unification into the 
CMSSM, is due to the generation of sizeable SUSY corrections to nib [HH which have 
the same sign as the MSSM parameter /x. Consequently, for /i > 0, the tree-level mfe(Mz), 
which is predicted (the sign of /x in this reference is opposite to the one adopted here) 
from Yukawa unification already above the 95% confidence level (c.l.) experimental range 
2.684 — 3.092 GeV, receives large positive corrections. This range is derived from the 
95% c.l. range for mb(mb) in the MS renormalization scheme which is found (see also 
Ref. 0) to be 3.95 - 4.55 GeV. This is evolved |[ up to Mz via the three-loop WS 
renormalization group (RG) equations with as{Mz) = 0.1185 and then converted 
to the DR scheme using the appropriate one-loop factor |1^. For /i < 0, the tree- 
level mb{Mz) is smaller and close to the upper edge of the above range. However, 
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the radiative corrections are now negative and drive mh{Mz) below this range. The 
discrepancy is, though, considerably smaller than in the /i > case. 

We see that, for both signs of /i, the hypothesis of exact Yukawa unification leads to an 
unacceptable 6-quark mass. However, we are not obliged to abandon Yukawa unification 
altogether. We can rather modestly correct it by including an extra SU{A)c non-singlet 
Higgs superfield with Yukawa couplings to the quarks and leptons. The Higgs SU{2)l 
doublets contained in this superfield can naturally develop [0 subdominant vacuum 
expectation values (vevs) and mix with the main electroweak doublets which are assumed 
to be SU{4:)c singlets and form a SU{2)r doublet. This mixing can, in general, violate 
SU{2)r. Consequently, the resulting electroweak Higgs doublets hf" , h'^ do not form 
a SU{2)r doublet and also break 5'f/(4)c. The required deviation from exact Yukawa 
unification is expected to be more pronounced in the /x > case. Despite this, we choose 
to study here this case since, for yU < 0, the present experimental data |jl2[ on the inclusive 
decay 6 — > 57 restrict the sparticle masses to considerably higher values and, thus, 
this case is phenomenologically less interesting. Moreover, the recent results |T^ on the 
muon anomalous magnetic moment also imply heavy sparticles for /i < (see below). 

We will construct here a concrete SUSY GUT model which naturally leads to a mod- 
est deviation from Yukawa unification allowing an acceptable mb(M^) even with uni- 
versal boundary conditions. (For models which violate universality rather than Yukawa 



unification see Ref. We will then show that this model possesses a wide range of pa- 

rameters which is consistent with all the phenomenological and cosmological constraints. 
We consider only the /i > case, which is experimentally more attractive. 

In Sec.|I|, we construct a SUSY GUT model which is based on the PS gauge group and 
provides, in a natural way, a suppressed violation of Yukawa unification. Variants of this 
model which can yield bigger deviations from Yukawa unification are also presented. In 



Sec.|T|, we concentrate on one of these variants which can violate Yukawa unification by 
an amount that is adequate for > 0. We then describe the resulting MSSM under the 
assumption of universal boundary conditions and introduce the various phenomenological 
and cosmological requirements which restrict its parameter space. In Sec.^, we study 
the range of parameters which is compatible with all these requirements in this particular 
CMSSM. Finally, in Sec.|V|, we summarize our conclusions. 
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II. THE SUSY GUT MODEL 



We take the SUSY GUT model of Ref. |T5[ (see also Ref. [|T6|) as our starting point. 
This is based on the PS gauge group Gps, which is the simplest gauge group that can 
lead to Yukawa unification. The 'matter' superfields are Fj = (4, 2, 1) and = (4, 1, 2) 
{i = 1,2,3), while the electroweak Higgs doublets belong to the superfield h = (1,2,2). 
So, all the requirements for exact Yukawa unification are fulfilled. The breaking of Gps 
down to the standard model (SM) gauge group Gs is achieved by the superheavy vevs 
(~ Mgut) of the right handed neutrino type components of a conjugate pair of Higgs 
superfields H'^ = (4, 1, 2), H'^ = (4, 1, 2). The model also contains a gauge singlet S which 



triggers the breaking of Gps, a SU{A)c 6-plet G = (6, 1, 1) which gives IjT^] masses to the 
right handed down quark type components of H^, H'^, and a pair of gauge singlets N, 
N for solving the fi problem of the MSSM via a Peccei-Quinn (PQ) symmetry. In 
addition to Gps, the model possesses two global U{1) symmetries, namely a PQ and a R 
symmetry, as well as a Z^^ symmetry ('matter parity') under which F, F'^ change sign. 
For details on the charge assignments, the full superpotential and the phenomenological 



and cosmological properties of this model, the reader is referred to Ref. [|T5 . 

A moderate violation of Yukawa unification can be accommodated in this model by 
adding a new Higgs superfield h' = (15, 2, 2) with Yukawa couplings FF'^h' . Actually, this 
is the only representation, besides (1,2,2), which possesses such couplings to the fermions. 
The existence of these couplings requires that the quantum numbers of h' coincide with 
the ones of h. So, its PQ and R charges are PQ = 1 and R = respectively. In order 
to give superheavy masses to the color non-singlet components of h', we need to include 
one more Higgs superfield h' = (15, 2, 2) with the superpotential coupling h'h', whose 
coefficient is of the order of Mgut- The field h' then has PQ = — 1 and R = 1. The full 
superpotential which is consistent with all the symmetries contains, in addition to the 



couplings mentioned above and in Ref. [|1^], the following extra terms too 



FFH^H^hh', FFH^H^h'h', FFH^H^hh', FFH^H^h'h', {H^fh'h, {Hyh'h, 

H'^H^h'h, N^H'^Yhh', N\H^yhh', N^H^H^hh', N^h'h'. (1) 

Note that all the superpotential terms can be multiplied by arbitrary powers of the 
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combinations {H^Y, {H'^)'^, H'^H'^. The first two couplings in Eq.(|I]) (as well as all 
the 'new' couplings containing {H'^Y, {H^Y) give rise to additional baryon and lepton 
number violation. However, their contribution to proton decay is subdominant to the 
one discussed in Ref. [|T5| and proton remains practically stable. 

Let us now discuss how the introduction of h', h' can lead to a moderate correction 
of exact Yukawa unification. The Higgs field h' contains two color singlet SU{2)l dou- 
blets h'2 which, after the breaking of Gps to Gs-, can mix with the corresponding 
doublets hi, h2 in h. This is mainly due to the terms h'h', H'^H^h'h. The latter, being 
non-renormalizable, is suppressed by the string scale Ms ~ 5 x 10^^ GeV. Actually, 
H^H^h'h corresponds to two independent couplings. One of them is between the SU{2)r 
singlets in H^W^ and h'h, and the other between the SU{2)Ft triplets in these combina- 
tions. The singlet (triplet) coupling, after the breaking of Gps to G5, generates SU{2)r 
invariant (violating) bilinear terms between the doublets in h' and h. So, we finally ob- 
tain two bilinear terms h'^hi and h'2h2 with different coefficients, which are suppressed by 
Mgut/Ms- These terms together with the terms h[h[ and h'2h'2 from h'h', which have 
equal coefficients, generate different mixings between hi, h'l and /i2, h'2. Consequently, the 
resulting electroweak doublets hl^, hf" contain SU{A)c violating components suppressed 
by Mgut/Ms and fail to form a SU{2)u doublet by an equally suppressed amount. So, 
Yukawa unification is moderately violated. Unfortunately, as it turns out, this violation 
is not adequate for correcting the 6-quark mass for /i > 0. 

In order to allow for a more sizable violation of Yukawa unification, we further extend 
the model by including a superfield = (15, 1, 3) with the coupling (j)h'h. The field is 
then neutral under both the PQ and R symmetries. To give superheavy masses to the 
color non-singlets in 0, we introduce one more superfield = (15, 1, 3) with the coupling 
00, whose coefficient is of order Mqut- The charges of are PQ = 0, i? = 1. Additional 
superpotential couplings which are allowed by the symmetries of the model are: 

^{H'Y, ^{H'Y, ^H'H", (PN^hh'. (2) 

Multiplications of the superpotential terms by 0, {H'^Y, {H^Y, H'^H^ are generally al- 
lowed with the exception of the multiplication of the Gps singlets S, N'^N'^, N'^h^ by a 
single 0, which is a S'f/(4)c 15-plet. 
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The terms 4>4> ^^^^ 4>H'^H'^ imply that, after the breaking of Gps to Gs, the field 
acquires a superheavy vev of order Mqut- The coupling (t)h'h then generates SU{2)r 
violating unsuppressed bilinear terms between the doublets in h! and h. These terms can 
certainly overshadow the corresponding ones from the non-renormalizable term H^H'^h'h. 
The resulting SU{2)f> violating mixing of the doublets in h and h' is then unsuppressed 
and we can obtain stronger violation of Yukawa unification. 

Instead of the pair 0, one can use a pair of SU{2)r singlet superfields (f)' = (15, 1, 1), 
(f)' — (15, 1, 1) with the same PQ and R charges and the couphngs (f)'h'h and (f)'(f)'. In 
this case, the mixing is SU{2)r invariant and the 'asymptotic' relation ht — hi, is not 
violated. The additional superpotential couplings allowed by the symmetries are: 



Note the absence of the couplings (f)'{H'^)^, (f)'{H'^y which are SU{2)r triplets. Multipli- 
cations of the superpotential terms by 0', {H^Y, {H'^) K H' H'' arc generally allowed with 
the exception of the multiplication of S, N'^N'^, N'^h'^ by a single 0', which is a SU{4:)c 
15-plet, or a single 4)'{H^Y, 4>'{H^Y, which are SU{2)r triplets. 

Finally, one could introduce both the pairs (j), (f) and 0', 0' at the same time. The 
allowed superpotential terms include all the above mentioned terms and 



Multiplications of the superpotential terms by 0, 0', {H'^Y, {H'^Y^ H'^H'^ are generally 
allowed with the exception of the multiplication of N'^N^, N'^h'^ by a single factor 
of 0, 0', 00', which arc 5'[/(4)e 15-plets, or a single factor of (t){H^Yi 4''{H'^Y- Also, 
multiplication of 0'0' by a single is not allowed. 

To further analyze the mixing of the doublets in h and h', we must first define properly 
the relevant couplings h'h', 4>h'h and (p'h'h. The superfield h — (1,2,2) is written as 



The color singlet components of the fields h', h' can be similarly represented. Under 

Uc G SU{4:)c, Ul & SU{2)l and Ur G SU{2)r, the relevant fields transforms as: 



^'{H'Y, ^'{H'Y, ^'H'H", <f>'N^hh'. 



(3) 



0'0(^^)^ 0'0(i^^)^ 0^. 



(4) 
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F ^ FUlil,, ^ UliliF^ h ^ ULhilR, h' ^ U^UJi'IIrUI h' ^ UlULh'UnU,, 

(f) ^ u,Ur<puIuI <p' u,<p'ul (6) 

where the transpose of a matrix is denoted by tilde. The Yukawa couphngs are FhF'^, 
Fh'F'^. From Eq.(^) and the identities Ul^Ul = Ur^Ur = e (e is the 2x2 antisymmetric 
matrix with ei2 = 1) which follow from det Ul = det Ur = 1, one finds that he URheUj^, 
h'e UJlJi'tU^Ul, h'e UMnh'tUlUl We see that tT{h'eh'e), iiiji'tfht), tT{h'e(f)'he) 
(the traces are taken with respect to the SU{A)c and SU{2)l indices) are invariant under 
Gps- They correspond to the 'symbolic' couplings h'h', (ph'h, (p'h'h respectively. 

After the breaking of Gps to Gs, the fields 0, 0' acquire vevs (0), (0') ~ Mqut- 
Substituting them by these vevs in the above couplings and using Eq.@, we obtain 



tT{h'eh'e) = h-^eh'^ + h\eh'^ + ■ ■ ■ , (7) 



id') - ~ id') ~ - 

tr(/i'e0'/ie) = ^ii{h'ehe) = ^{h^eh2 + hieh'^), (9) 
V 2 V 2 

where only the colorless components of h' and h' are shown in the right hand side of 
Eq.(0) and = diag(l, —1). The bilinear terms between hi, h[, h[ and /i2, h2, h'2 which 
appear in Eqs.(0), (|) and (H) turn out to be the dominant bilinear terms between these 
doublets. Collecting them together, we obtain 

mh^e{h'2 + 0:2^2) + rnih'^ + aihi)eh'2i (10) 

where m is the superheavy mass parameter which multiplies the term in Eq.(|^) and 
«! = (— p(0)+g(0'))/-\/2m, a2 = (p(0)+g(0'))/v^m with p and g being the dimensionless 
coupling constants which correspond to the SU{2)r triplet and singlet terms in Eqs.(H) 
and (ID respectively. Note that ai, 0:2 are in general complex. So, we get two pairs of 
superheavy doublets with mass m. They are predominantly given by 

K , and , h'2. (11) 

(l + |a2P)5 (l + l"iP)^ 

The orthogonal combinations of hi, h[ and /12, h2 constitute the electroweak doublets: 

hr = and hT = . (12) 

(l + |ai|2)l (l + lttsH^ 



We see that, although hi, h2 and h[, h'2 form SU{2)r doublets, this is, in general, not 
true for , since ai, 02 can be different. However, if we remove from the model 
the SU{2)r triplets 0, (and, thus, the SU{2)r triplet coupling in Eq.(|])), we obtain 
«! = a2 and hl^, h*^ do form a SU{2)ji doublet. The minimal model which can provide 
us with an adequate violation of t — 6 Yukawa unification for /i > is the one including 
0, but not the 5'f/(2)/j singlets 0', 0'. This model, which we choose to study in detail 
here, yields ai = —a^- Inclusion of both 0, and 0', 0' can lead to an arbitrary relation 
between a\ and 02 and, thus, introduces an extra complex parameter. Finally, due to 
the presence of h!^, h^, the electroweak doublets are generally non-singlets under SU{'i)c 
and h — T Yukawa unification can also be violated. 

The doublets in Eq.(ll) must have zero vevs, which implies that {h'^) = —ai{hi), 
{h'2) = -a2{h2). Eq.(|12D then gives (V) = (1 + \ai\^y/^{h), {h^^) = (1 + |«2p)'/'(/i2). 
From the third generation Yukawa couplings y^sF^hF^, 2y'^^F3h' F^ , we obtain 

1 



mt 



|2/33(/i2) +2/33(^2)1 = l(?/33 - 2/33"2)(/i2)| 



-y 33 [1^2 I 



(1 + |a2 



2U ' 



(13) 



where p = ^33/2/33 and can be taken positive by appropriately readjusting the phases 
of h, h' . Note that, in the SU{A)c space, the doublets in h' are proportional to 
diag(l/2\/3, l/2v^, l/2v^, — v^/2), which is normalized so that the trace of its square 
equals unity. Thus, to make 7/33 directly comparable to 2/33, we included a factor of two 
in defining the corresponding Yukawa coupling. The masses mj,, irir are, similarly, found: 



1 + v^pai 



'1 + kil 



-y33{hT) 



(14) 



From Eqs. (p!3D and (|lj), we see that the exact equality of the 'asymptotic' Yukawa 
couplings ht, hh, h^. is now replaced by the quasi- unification condition: 



hi : h: K 



pa2 



/V3 



[I +1^2 



|2U 



1 — poi/a/S 



(1 + |ai 



1 + \^pai 



(l + |ai|2)- 



(15) 



This condition depends on two complex {ai, 02) and one real (p > 0) parameter. 

Note that the mixing between hi, h[ (i.e., ai 7^ 0) and the fact that the Higgs 
superfield h' possesses Yukawa couplings to the 'matter' superfields (i.e., p 7^ 0) are 
crucial for violating b — r Yukawa unification, which is though not affected by the mixing 
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between /i2, h'2 (i.e., the value of a2)- On the contrary, violation of t— 6 Yukawa unification 
requires that ai 7^ 0:2, which can be achieved only in the presence of a SU{2)r triplet 
bilinear term between h' and h. In summary, the minimal requirements for full violation 
of Yukawa unification are p, ai 7^ 0, ai 7^ a2- In the minimal model (with 0, but not 



b') which we will study here, ai = —a2 and, thus, Eq.(|l5D takes the simple form 



/li : /ife : /i^ = |1 + c| : |1 - c| : |1 + 3c|, (16) 

where c = pai/-\/3. This 'asymptotic' relation depends on a single complex parameter. 
For simplicity, we will restrict our analysis to real values of c only. 

For completeness, we also give the 'asymptotic' relation between the Yukawa couplings 
in the model without 0, 0, 0', 0', which, although not suitable for /i > 0, may be 
adequate for p < 0. Under Gps-, the superfields H'^, H'^ transform as: H'^ H'^UrUc, 
tJltjiH^. Thus, HH^ UcUrIi'' H^U^Ul and ti{h'eHH^he) is invariant under 
Gps containing both the SU{2)p singlet and triplet 'symbolic' couplings H'^H'^h'h. In the 
SU{2)r space, {H = diag(wg, 0) = vlao/2 + vla^/2, where vq = Mgut/9gut, with 

dGUT being the GUT gauge coupling constant, and ctq is the unit 2x2 matrix. Moreover, 
in the S'f/(4)c space, {H ){H^) = diag(0, 0, 0, fg). So, replacing H^, H'^ by their vevs in 
the singlet and triplet couplings in tr{h'eH H'^he) (with dimensionless coefficients equal 
to q' and p'), we obtain {—\/?)vl/ A)ii{h'ehe) and {—\/?)VQ/A)ii{h'ea'ihe). We, thus, again 
end up with Eq.(|l|) but with ai = -^/3v^{q' - p') /AmMs, 02 = -V^v'^{q' + p')/4mMs, 
which are suppressed by Mqut / Ms- 



Ill. THE RESULTING MSSM 

We will now concentrate on the minimal model which includes 0, but not 0', 0'. This 
model, below Mqut, reduces to the MSSM supplemented by the 'asymptotic' Yukawa 
coupling quasi-unification condition in Eq.(|l6D, where c is taken real for simplicity. We 
will assume universal soft SUSY breaking terms at Mqut, i-e., a common mass for all 
scalar fields mo, a common gaugino mass M1/2 and a common trilinear scalar coupling 
Aq. So the resulting MSSM is actually the CMSSM. Furthermore, we will concentrate 
on the fi > case for reasons which we already explained. 
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We will closely follow the notation as well as the RG and radiative electroweak break- 
ing analysis of Ref. [jT^ for the CMSSM with the improvements of Ref. (recall that 
the sign of /i in these references is opposite to the one adopted here). These improve- 
ments include the employment of the full one-loop corrections to the effective potential 
for the electroweak breaking and to certain particle masses taken from Ref. |^|. They 
also include the incorporation of the two-loop corrections to the CP-even neutral Higgs 
boson mass matrix by using FeynHiggsFast and the introduction [^] of a variable 
common SUSY threshold Msusy = {f^ii'^hY^'^ (^1,2 are the stop mass eigenstates) where 
the effective potential is minimized, itla (the CP-odd Higgs boson mass) and /i are eval- 
uated, and the MSSM RG equations are replaced by the SM ones. Note that, at Msusy, 



the size of the one- loop corrections to and is [22| minimal and, thus, the accuracy 
in the determination of these quantities is maximal. 

For any given mb{Mz) in its 95% c.l. range (2.684 — 3.092 GeV), we can determine the 
parameters c and tan (3 at Msusy so that the 'asymptotic' condition in Eq.(|16|) is satisfied. 
We use fixed values for the running top quark mass mt{mt) = 166 GeV and the running 
tau lepton mass mr{Mz) = 1.746 GeV and incorporate not only the SUSY correction 
to the bottom quark mass but also the SUSY threshold correction to mr{MsusY) from 
the approximate formula of Ref. pn|. This correction arises mainly from chargino/tau 
sneutrino (ur) loops and, for /i > 0, leads p to a small reduction of tan/3. 

After imposing the conditions of gauge coupling unification, successful electroweak 
breaking and Yukawa quasi- unification in Eq.(|l^), we are left with three free input pa- 
rameters mo, Ml/2 and Aq (see Refs. [|,0]). In order to make the notation physically 
more transparent, we replace mo and M1/2 equivalently by the mass musp (or ^x) '^^ 
the lightest supersymmetric particle (LSP), which turns out to be the lightest neutralino 
(x), and the relative mass splitting Af^ = {nif^ —m^)/m^ between the lightest stau mass 
eigenstate (f2) and the LSP. We will study the parameter space of this model which is 
compatible with all the available phenomenological and cosmological constraints. 

An important constraint is obtained by considering the inclusive branching ratio 
BR(6 —>■ 57) of the decay process b —>■ 57. The present best average of this branching 
ratio, which is found from the available experimental data []T^|, is 3.24 x 10~^ with an 
experimental error of ±0.38 x 10~^ and an asymmetric 'theoretical' error due to model 
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dependence of [+0.26, —0.23] x 10~^ (assuming no correlation between the experimental 
systematics, which should not be too far from reality). 

We calculate BR(6 57) using the formalism of Ref. where the SM contribu- 
tion is factorized out. This contribution includes the next-to-leading order (NLO) QCD 
and the leading order (LO) QED corrections. The charged Higgs boson contribution to 



BR(6 — > 57) is evaluated by including the NLO QCD corrections from Ref. [^. The 
dominant SUSY contribution includes the NLO QCD corrections from Ref. [^], which 
hold for large tan/? (see also Ref. [^). The error in the SM contribution to BR(6 37) 
is made of two components. The first is due to the propagation of the experimental 



errors in the input parameters and is about ±0.23 x 10 . The second originates 



from the dependence on the renormalization and matching scales and is [28| about ±6% 



if the NLO QCD corrections are included. These errors remain p8[ basically the same 
even after including the charged Higgs and SUSY contributions corrected at NLO. 

In order to construct the 95% c.l. range of BR(6 — > S7), we first add in quadrature 
the experimental error (±0.38 x 10^^) and the error originating from the input param- 
eters (±0.23 X 10"'^). This yields the overall standard deviation (a). The asymmetric 
error coming from model dependence ([±0.26,-0.23] x 10^"^) and the error from scale 
dependence (6%) are then added linearly [^ on both ends of the ±2 — cr range. For 
simplicity, we take a constant error from scale dependence, which we evaluate at the 
central experimental value of BR(6 57). The 95% c.l. range of this branching ratio 
then turns out to be about (1.9 — 4.6) x 10^^. 

The recent measurement [|l^] of the anomalous magnetic moment of the muon = 
{qh — 2)/2 provides an additional significant constraint. The deviation of from its 
predicted value in the SM, Sa^, is found to lie, at 95% c.l., in the range from —6 x 10^^*^ 
to 58 X lO"^''. This range is derived using the calculations (see e.g., Ref. |^0[) of in the 
SM which are based on the evaluation of the hadronic vacuum polarization contribution 



of Ref. However, we take here the recently corrected [^ sign of the pseudoscalar 
pole contribution to the light-by- light scattering correction to a^. This corrected sign 
reduces considerably the discrepancy between the SM and the measured value of a^. It 
also relaxes the restrictions on the parameter space of the CMSSM. In particular, the 
H < case, which leads to negative (5a^, is no longer disfavored. However, the sparticles, 
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in this case, cannot be as light as in the jj, > case, where Sa^ > 0. The calculation of 
Sa^ in the CMSSM is performed here by using the analysis of Ref. PB . 



Another constraint results from the requirement that the relic abundance ^Ilsp h"^ of 
the LSP in the universe does not exceed the upper limit on the cold dark matter (CDM) 
abundance which is derived from observations {^Ilsp is the present energy density of the 
LSP over the critical energy density of the universe and h is the present value of the 
Hubble parameter in units of 100 km sec~^ Mpc~^). From the recent results of DASI 
3^ , one finds that the 95% c.l. range of flcDM h"^ is 0.06 — 0.22. Therefore, we require 
that Qlsp does not exceed 0.22. 

Here, the LSP (x) is an almost pure bino. Its relic abundance will be calculated by 



micrOMEGAs |35|, which is the most complete code available. (A similar calculation has 
appeared in Ref. [0.) It includes all the coannihilations [^] of neutralinos, charginos, 
sleptons, squarks and gluinos. The exact tree-level cross sections are used and are accu- 
rately thermally averaged. Also, poles and thresholds are properly handled and one-loop 
QCD corrected Higgs decay widths |33| are used, which is the main improvement pro- 



vided by Ref. [Q. The SUSY corrections to these widths are, however, not included. 
Fortunately, in our case, their effect is much smaller than that of the QCD corrections. 

In order to have an independent check of micrOMEGAs, we also use the following al- 
ternative method for calculating Qlsp in our model. In most of the parameter space 
where coannihilations are unimportant, ^Ilsp h"^ can be calculated by using DarkSUSY 



40|| . This code employs the complete tree- level (and in some cases one-loop) cross sec- 



tions of the relevant processes with properly treating all resonances and thresholds and 
performs accurate numerical integration of the Boltzmann equation. Its neutralino anni- 
hilation part is in excellent numerical agreement with the recent exact analytic calculation 



of Ref. I^!]]. Moreover, it includes all coannihilations between neutralinos and charginos 
which are, however, insignificant for our model. Its main defect is that it uses the tree- 
level Higgs decay widths. This can be approximately corrected if, in evaluating the Higgs 
decay widths, we replace mb{mb) by nib at the mass of the appropriate Higgs boson in 



the couplings of the 6-quark to the Higgs bosons (see Ref. p5[]). 

In the region of the parameter space where coannihilations come into play, the next- 
to-lightest supersymmetric particle (NLSP) turns out to be the f2 and the only relevant 
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coannihilations are the bino-stau ones []T9|,0. In this region, which is given by A^^ < 0.25, 
we calculate flisp h"^ by using an improved version of the analysis of Ref. (the sign 
of /i in this reference is opposite to the one adopted here). This analysis, which has 
been applied in Refs. [§,0], includes bino-stau coannihilations for all tan/3's (see also 
Refs. 0,0) and is based on a series expansion of the thermally averaged cross sections 
in = Tp/m^, with Tp being the freeze-out temperature. This expansion is, however, 
inadequate when we have to treat resonances, and these are crucial for bino annihilation. 



So, we need to improve the bino annihilation part in Ref. |T9[. This can be achieved by 
evaluating the corresponding expansion coefficients a^^ and h^^ not as in this reference 
but as follows. From DarkSUSY, we find the values of Vty. h'^ and Xp which correspond to 
essentially just bino annihilation. Using the appropriate formulas of Ref. [1^, we then 
extract the 'improved' expansion coefficients a^.^. and b^^. 

There is one more improvement which we need to do in Ref. ||19[ (and was already 
used in Ref. [^) to make it applicable to the present case. The cross sections of the 
processes f2f2 W~^W~ , H^H^ were calculated without including the tree-graphs 
with Vr exchange in the t-channel |^6[. Thus, the contribution to a^jf^ from the process 
T'2T'2 ~^ W^W~ appearing in Table II of Ref. ||T^ should be corrected by adding to it 



,2 yl - m^)"m4 

i2 I ^2 



(17) 



/o Q - 2 , ^ 4. \ I 9HgHW+W- , 9hghW+W- , 2 

-(2 - 3m^^ + m^,) — \ — + gf^f^w+w-m^. 

with gf^i^^w± = gs9/V2. Also, the contribution to af^f* from f2r2 H^H^ is now given 
by Eq.(26) of Ref. [l^ with the expression in the last parenthesis corrected by adding to 
it the quantity gf^u-rH^/i^'iii ~ " -'-)' where 



gf2VrH± = fl'[("^vy — m^tan/?)se + (A^-tan/? -|- jj)'mrCe\/ y2m 

(with the present sign convention for /i). The resulting correction to VLy. h'^, in the case 
of the process 7^2X2 W^W~ , varies from about 8% to about 1% as Af^ increases from 
to 0.25. On the contrary, the correction from f2r| H^H~ is negligible. 

We find that the alternative method for calculating the neutralino relic abundance in 
our model, which we have just described, yields results which are in excellent agreement 
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with micrOMEGAs. In practice, however, we use the code micrOMEGAs since it has the 
extra advantage of being much faster. 

We will also impose the 95% c.l. LEP bound on the lightest CP-even neutral Higgs 
boson mass ruh > 114.1 GeV |^^. In the CMSSM, this bound holds almost always for all 
tan/5's, at least as long as CP is conserved. Finally, for the values of tan/3 which appear 
here (about 60), the CDF results yield the 95% c.l. bound uia > HO GeV p . 



IV. THE ALLOWED PARAMETER SPACE 



We now proceed to the derivation of the restrictions which are imposed by the various 
phenomenological and cosmological constraints presented in Sec.|T| on the parameters of 
the CMSSM (with ^ > 0) supplemented with the Yukawa quasi-unification condition in 
Eq.(|16[) (with < c < 1). The restrictions on the rriLSP — plane, for = and 
with the central value of as{Mz) = 0.1185, are shown in Fig.l. The dashed (dotted) 
lines correspond to the 95% c.l. lower (upper) experimental bound on rrih^Mz) which 
is 2.684 GeV (3.092 GeV) (see SecD, while the solid lines correspond to the central 
experimental value of mb{Mz) = 2.888 GeV. We will follow this convention in the 
subsequent four figures too. From left to right, the dashed (dotted) hues depict the 95% 
c.l. lower bounds on misp from the constraints > 110 GeV, BR(6 — > 37) > 1.9 x 10~^ 
and dttfj, < 58 x 10^^°, and the 95% c.l. upper bound on rriLSP from VLlsp h"^ < 0.22. 
The constraints BR(6 57) < 4.6 x 10~^ and > —6 x 10^^" do not restrict the 
parameters since they are always satisfied for ^ > 0. The left solid line depicts the lower 
bound on itllsp from rrih > 114.1 GeV which does not depend much on mi){Mz), while 
the right solid line corresponds to Qlsp h"^ = 0.22 for the central value of mb{Mz)- 

We observe that the lower bounds on ttilsp are generally not so sensitive to the 
variations of mb{Mz) within its 95% c.l. range. The most sensitive of these bounds 
is the one from > 110 GeV, while the bound from rrih > 114.1 GeV is practically 
mfe(M^)-independent. Actually, this bound overshadows all the other lower bounds for all 
Afj's and, being essentially A^j-independent too, sets, for Aq = and as{Mz) = 0.1185, 
an overall constant 95% c.l. lower bound on rriLSP of about 138 GeV. Contrary to the 
lower bounds, the line from the LSP relic abundance is extremely sensitive to the value 
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of mf,{Mz)- In particular, its almost vertical part is considerably displaced to higher 
rriLSp'^ as mi,{Mz) decreases. We will explain this behavior later. 

In Fig. 2, we depict and Msusy versus rriLSP for various A^^j's, = and the 
central values of mb{Mz) and as{Mz)- We see that is always smaller than 2mLSP 
but close to it. Thus, generally, the neutralino annihilation via the s-channel exchange 
of an A-boson is by far the dominant (co) annihilation process. We also observe that, as 
fTT-LSP or increase, we move away from the A-po\e, which thus becomes less efficient. 
As a consequence, Qlsp increases with rriLSP or /S.f^ (see Fig. 3). 

In Fig. 3, we show VLlsp as a function of misp for various A^j's, Aq = 0, nibiMz) = 
2.888 GeV and asiMz) = 0.1185. The sohd lines are obtained by using micrOMEGAs. The 



crosses are from the alternative method described in Sec.|lTl|, which combines DarkSUSY 
with the bino-stau coannihilation calculation of Ref. In order to mimic, in DarkSUSY, 
the effect of the one-loop QCD corrections to the Higgs decay widths, we can generally 
use the 6-quark mass (including the SUSY corrections) evaluated at rriA since the A-pole 
is by far the most important. As it turns out, mhimA) depends very weakly on the 
values of the input parameters rriLSP and Af^, and thus the uniform use of a constant 
mean value of mi,{'mA) should be adequate. We find that excellent agreement between 
micrOMEGAs and our alternative method is achieve if, in DarkSUSY, we use a constant 
default value for the 6-quark mass which is equal to about 2.5 GeV (see Fig. 3). 

The importance for our calculation of the one-loop QCD corrections to the Higgs 
decay widths can be easily concluded from Fig. 4, where we show VLlsp h'^ versus rriLSP for 
Af2 = I, Aq = and the central values of mb{Mz), as{Mz)- The LSP relic abundance is 
calculated by micrOMEGAs with (thick solid line) or without (faint solid line) the inclusion 
of the one-loop QCD corrections to the Higgs widths. We see that the results differ 
appreciably. For comparison, we also draw ^Ilsp h"^ calculated by our alternative method 
with (thick crosses) or without (faint crosses) the above 'correction' of the 6-quark mass. 
The agreement with micrOMEGAs is really impressive in both cases. 

In Fig. 5, we present rriA and Msusy as functions of ttilsp for A^^ = 1, Aq = and 
with the lower, upper and central values of mb{Mz) for as{Mz) = 0.1185. We see that 
tua increases and approaches 2mLSP as mb{Mz) decreases. As a consequence, the effect 
of the A-pole is considerably enhanced and ^Ilsp is reduced. This explains the large 
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displacement of the vertical part of the Qlsp = 0.22 line in Fig.l at low mft(M^)'s. 

As becomes small and approaches zero, bino-stau coannihilations become very 
important and strongly dominate over the pole annihilation of neutralinos. This leads 
to a very pronounced reduction of the neutralino relic abundance and, consequently, to 
a drastic increase of the upper bound on itllsp- We, thus, obtain the almost horizontal 
parts of the ^Ilsp h"^ = 0.22 lines in Fig.l. 

The shaded area depicted in Fig.l, which lies between the = 114.1 GeV line 
and the Qlsp = 0.22 line corresponding to the central value of mfe(M^), constitutes 
the allowed range in the ttilsp — plane for Aq = 0, mb{Mz) = 2.888 GeV and 
as{Mz) = 0.1185. For practical reasons, we limited our investigation to Af^ < 2 and 
^LSP < 500 GeV. We see that allowing mi,{Mz) to vary within its 95% c.l. experimental 
range considerably enlarges the allowed area in Fig.l. On the other hand, fixing mb{Mz) 
to its upper bound yields the smallest possible allowed area for Aq = and as{Mz) = 
0.1185, which is, however, still quite wide. In any case, we conclude that, for Aq = 0, 
there exists a natural and wide parameter space which is consistent with all the available 
phenomenological and cosmological requirements. 

In the allowed (shaded) area of Fig.l which corresponds to the central value of 
mb{Mz), the parameter c (tan/3) varies between about 0.15 and 0.20 (58 and 59). For 
mb{Mz) fixed to its lower or upper bound, we find that, in the corresponding allowed 
area, the parameter c (tan/5) ranges between about 0.17 and 0.23 (59 and 61) or 0.13 and 
0.17 (56 and 58). We observe that, as we increase mb{Mz), the parameter c decreases 
and we get closer to exact Yukawa unification. This behavior is certainly consistent with 
the fact that the value of mh{Mz) which corresponds to exact Yukawa unification lies well 
above its 95% c.l. range. The LSP mass is restricted to be higher than about 138 GeV 
for y4o = and as{Mz) = 0.1185, with the minimum being practically -independent. 
At this minimum, c ^ 0.16 - 0.20 (c ^ 0.13 - 0.23) and tan/5 ^ 59 (tan/? ^ 58 - 61) for 
mb{Mz) = 2.888 GeV (mfe(Mz) = 2.684 - 3.092 GeV). 

The 'asymptotic' Yukawa quasi-unification condition in Eq.(0) with < c < 1 yields 

su h-ht hr-ht 2c 

bh= ; = — = - (18) 

ht ht l + c 

This means that the bottom and tau Yukawa couplings split from the top Yukawa cou- 
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pling by the same amount but in opposite directions, with hi, becoming smaller than 
ht- For Aq — and with the central values of mb{Mz) and as{Mz), the sphtting 6h 
ranges from about 0.26 to 0.33. Allowing mi,{Mz) to vary in its 95% c.l. range, how- 
ever, we obtain a larger range for 5h which is about 0.22 — 0.38. At the minimum of 
niLSP ~ 138 GcV , 5h ^ 0.28 - 0.33 {dh ^ 0.23 - 0.38) for rribiMz) = 2.888 GeV 
{mb{Mz) = 2.684 - 3.092 GcV). 

For simplicity, as{Mz) was fixed to its central experimental value (equal to 0.1185) 
throughout our calculation. We could, however, let it vary in its 95% c.l. experimental 
range 0.1145 — 0.1225, which would significantly widen the 95% c.l. range of mb{Mz) 
to 2.616 — 3.157 GeV, with the lower (upper) bound corresponding to the upper (lower) 
bound on as{Mz). This would lead to a sizable enlargement of the allowed area due to 
the sensitivity of the LSP relic abundance to the 6-quark mass. 

In Figs. 6, 7 and 8, we present the restrictions on the uilsp — ^0/^1/2 plane for 
mb{Mz) = 2.888 GeV, a,(Mz) = 0.1185, and A^^ = 0, 1 and 2 respectively We use 
solid, dashed, dot-dashed and dotted lines to depict the lower bounds on itllsp from 
niA > 110 GeV, BR(6 ^ 57) > 1.9 x 10~*, 5a^ < 58 x 10-^° and mh > 114.1 GcV 
respectively. The upper bounds on rriLSP from VIlsp < 0.22 are represented by double 
dot-dashed lines. Note that, in Fig. 6, the upper bound on itilsp does not appear since it 
hes at rriLSP > 500 GeV. This is due to the fact that, since the bino and stau masses are 
degenerate in this case, the bino-stau coannihilation is considerably enhanced leading to 
a strong reduction of the LSP relic abundance. In order to illustrate the quick reduction 
of ^LSP as approaches zero, we display in Fig.6 the bounds from ^lsp h"^ < 0.22 
(double dot-dashed hues) for Af^ = 0.1 and 0.03 too. The reduction of VIlsp h"^ by 
bino-stau coannihilation is also the reason for the fact that the allowed (shaded) areas in 
Figs. 6, 7 and 8 are more sizable for smaller A^^'s. The lower bound on itllsp overshadows 
all the other lower bounds in all cases and is only very mildly dependent on Afj. 

In Figs. 9 and 10, we present the restrictions on the uilsp — Afj plane for mb{Mz) — 
2.888 GeV, q;,(Mz) = 0.1185, and = 3 and - 3 respectively. We follow the 

same notation for the various bounds on rriLSP and the allowed areas as in Figs. 6, 7 and 
8. We observe that, for Aq/Mii2 — 3 (^o/-^i/2 = —3), the lower bound on rriLSP from 
the constraint BR(6 — > 57) > 1.9 x 10~^ is overshadowed by (overshadows) the bounds 
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from the requirements > 110 GeV and 6a^ < 58 x 10~^°. In all cases, however, the 
overall lower bound on misp is set by rrih > 114.1 GeV. 

For presentation purposes, we restricted our analysis to rriLSP < 500 GeV, < 2, 
—3 < Ao/Mi/2 < 3. The allowed ranges in Figs.l, 6, 7, 8, 9 and 10 constitute sections or 
boundaries of the part of the overall allowed parameter space for mi,{Mz) = 2.888 GeV 
and as{Mz) = 0.1185 which is contained in the investigated range of parameters. We 
found that, within this part of the overall allowed parameter space, tan/3 ranges between 
about 58 and 61. Also, the parameter c ranges between about 0.15 and 0.21 and, thus, the 
'asymptotic' splitting between the bottom (or tau) and the top Yukawa couplings varies 
in the range 26 — 35%. Moreover, the smallest value of rriLSP is about 107 GeV and is 
achieved at A0/M1/2 = —3 and practically any Af^. The corresponding 6h ^ 0.30 — 0.35 
and tan/5 60 — 61. Allowing mb{Mz) and as{Mz) to vary within their 95% c.l. ranges 
will further widen the overall allowed parameter space within the investigated range of 
parameters as well as the allowed ranges of the parameters c (or 6h) and tan/5. The 
minimal value of rriLSP is, though, not much affected by varying m[,{Mz)- Variations of 
as{Mz)-, however, lead to small (~ 10 GeV) fluctuations of the minimal rriLSP- 

We see that the required deviation from Yukawa unification, for yU > 0, is not so 
small. In spite of this, the restrictions from Yukawa unification are not completely lost 
but only somewhat weakened. In particular, tan j3 remains large and close to 60. Actually, 
our model is much closer to Yukawa unification than generic models where the Yukawa 
couplings can differ even by orders of magnitude. Also, the deviation from Yukawa 
unification is generated here in a natural, systematic, controlled and well-motivated way. 
Finally, recall that the required size of the violation of Yukawa unification forced us to 
extend the initial model of Sec.|| which could only provide a suppressed violation. 

V. CONCLUSIONS 

We constructed a class of concrete SUSY GUTs based on the PS gauge group which 
naturally lead to a moderate violation of 'asymptotic' Yukawa unification so that the 
6-quark mass can take acceptable values even with universal boundary conditions. For 
/X < 0, a suppressed deviation from Yukawa unification may be adequate, while, for 
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// > 0, a more sizable deviation is required. In the // < case, however, the sparticles 
are considerably heavier due to the constraints from 6 — > 57 and the muon anomalous 
magnetic moment. So, the /i > case is more attractive for experimenters. 

We considered a particular SUSY GUT from the above class with a deviation from 
Yukawa unification which is adequate for /i > 0. We then discussed the resulting MSSM 
under the assumption of universal boundary conditions and the various phenomenological 
and cosmological requirements which restrict its parameter space. They originate from 
the data on the inclusive branching ratio of 6 — >^ sj, the muon anomalous magnetic 
moment, the CDM abundance in the universe, and the masses and itia- 

The calculation of BR(6 — > 57) incorporates all the LO QED and NLO QCD cor- 
rections which hold for large values of tan/3, while the LSP contribution to flcDM is 
evaluated by using the code micrOMEGAs which includes all possible coannihilation pro- 
cesses, treats poles properly and uses the one-loop QCD corrected Higgs decay widths. 
We also employed an alternative method for estimating the LSP relic abundance and 
found excellent agreement with micrOMEGAs. 

We showed that, in the particular model with Yukawa quasi-unification considered, 
there exists a wide and natural range of CMSSM parameters which is consistent with 
all the above constraints. We found that, within the investigated part of the overall 
allowed parameter space, the parameter tan (5 ranges between about 58 and 61 and the 
'asymptotic' splitting between the bottom (or tau) and the top Yukawa couplings varies 
in the range 26 — 35% for central values of mb{Mz) and as{Mz). Also, the LSP mass 
can be as low as about 107 GeV. 
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FIGURES 




FIG. 1. Restrictions on the rriLSP — plane for Aq — 0, as{Mz) = 0.1185. From left to right, the 
dashed (dotted) lines depict the lower bounds on rriLSP from ttia > 110 GeV, BR(6 57) > 1.9 x 10~^ 
and 5a^ < 58 x 10"-^°, and the upper bound on niLSP from Qlsp h"^ < 0.22 for rUbiMz) = 2.684 GeV 
(3.092 GeV). The left (right) solid line depicts the lower (upper) bound on tulsp from ruh > 114.1 GeV 
{flLSP h? < 0.22) for mb{Mz) = 2.888 GeV. The allowed area for mb{Mz) = 2.888 GeV is shaded. 




FIG. 2. The mass parameters niA and Msusy as ftmctions of m,LSP for various values of Af^, which 
are indicated on the curves, and with Aq = 0, mb{Mz) = 2.888 GeV, ag{Mz) = 0.1185. 




FIG. 3. The LSP relic abmidancc Ulsp versus ttilsp for various Af2's (indicated on the curves) 
and with Aq = 0, mt{Mz) = 2.888 GcV, as{Mz) = 0.1185. The solid lines (crosses) are obtained by 
micrOMEGAs (our alternative method). The upper bound on ^Ilsp (=0.22) is also depicted. 




mj^p (GeV) 

FIG. 4. The LSP relic abundance ^lsp h'^ versus ttilsp for Af^ = I, Aq = 0, mb{Mz) = 2.888 GeV, 
as{Mz) = 0.1185. The thick (faint) solid line is obtained by micrOMEGAs with (without) the one-loop 
QCD corrections to the Higgs widths, while the thick (faint) crosses by our alternative method with 
(without) the 'correction' of the 6-quark mass. The upper bound on ^lsp h'^ (=0.22) is also depicted. 
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FIG. 5. The mass parameters tha and Msusy versus rriLSP for Af^ = 1, 
and with mb{Mz) = 2.684 GeV (dashed hnes), 3.092 GeV (dotted Unes) or 2.88^ 



= 0, a,(Mz) = 0.1185 
GeV (solid hnes). 
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FIG. 6. Restrictions on the rriLSP - MIM^ii plane for Af^ = 0, mb{Mz) = 2.888 GeV, as{Mz) = 
0.1185. The solid, dashed, dot-dashed and dotted lines correspond to the lower bounds on rriLSP from 
ruA > llOGeV, BR(6^ S7) > 1.9x10"'', 6a^ < 58x10"!° andm/i > 114.1 GeV respectively. The upper 
bound on rriLSP from flLSP < 0.22 does not appear in the figure since it lies at ttilsp > 500 GeV. 
The allowed area is shaded. For comparison, we also display the bounds from flLSP < 0.22 (double 
dot-dashed lines) for Af^ =0.1 and 0.03, as indicated. 
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FIG. 7. Restrictions on the rriLSP — ^0/^^1/2 plane for Af2 = 1, mb{Mz) = 2.888 GeV and 
as{Mz) = 0.1185. We use the same notation for the hnes which correspond to the various bounds on 
iTiLSP and for the aUowed area as in Fig. 6. The upper bound onniLSP from the cosmological constraint 
^LSP < 0.22 (double dot-dashed Une) now appears inside the figure. 
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FIG. 8. Restrictions on the ttilsp - ^0/^1/2 plane for Af^ = 2, mb{Mz) = 2.888 GeV and 
as{Mz) = 0.1185. We use the same notation for the hnes which correspond to the various bounds on 
iTiLSP and for the allowed area as in Fig. 7. 
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FIG. 9. Restrictions on the ttilsp - Af^ plane for An/Mi/2 = 3, mb{Mz) = 2.888 GcV and 
as{Mz) = 0.1185. We use the same notation for the hnes which correspond to the various bounds on 
mLSP and for the allowed area as in Fig. 7. 
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FIG. 10. Restrictions on the m.LSP ~ plane for A0/M1/2 = -3, mb{Mz) = 2.888 GeV and 
as{Mz) = 0.1185. We use the same notation for the lines which correspond to the various bounds on 
rriLSP and for the allowed area as in Fig. 7. The minimal value of thlsp in the investigated overall 
allowed parameter space is about 107 GeV for central values of mb{Mz) and as{Mz). As seen from this 
figure, this value is achieved at A0/M1/2 = —3 and practically any Af^ (see also Figs. 6, 7 and 8). 



